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Abstract. 

We study the real-time evolution of a charge-excited state in a one-dimensional e g -orbital 
degenerate Hubbard model, by a time-dependent density-matrix renormalization group method. 
Considering a chain along the z direction, electrons hop between adjacent 3z 2 —r 2 orbitals, while 
x 2 —y 2 orbitals are localized. For the charge-excited state, a holon-doublon pair is introduced 
into the ground state at quarter filling. At initial time, there is no electron in a holon site, while 
a pair of electrons occupies 3z 2 — r 2 orbital in a doublon site. As the time evolves, the holon 
motion is governed by the nearest-neighbor hopping, but the electron pair can transfer between 
3z 2 —r 2 orbital and x 2 —y 2 orbital through the pair hopping in addition to the nearest-neighbor 
hopping. Thus holon and doublon propagate at different speed due to the pair hopping that is 
characteristic of multi-orbital systems. 



1. Introduction 

Competing interaction among charge, spin, and orbital degrees of freedom plays a crucial role 
in the emergence of various types of quantum phases in strongly correlated electron systems 
with active orbital degree of freedom [H[2j[3]. The change of conditions, such as temperature, 
magnetic field, pressure, and chemical doping, gives rise to a dramatic change of transport and 
magnetic properties through a phase transition because of a subtle balance between multiple 
phases. 

Recent developments in femto-second laser technology have made it possible to investigate the 
dynamical control of the many-body electron state on an ultrafast time scale, the so-called photo- 
induced phase transition [3]. In general, the photo-irradiation can excite the system to a novel 
non-equilibrium state that cannot be accessed by simply changing temperature, since photon 
energy is much higher than thermal energy. A photo-excited non-equilibrium state eventually 
relaxes to the original equilibrium state. However, the system exhibits a transient behavior and 
sometimes it decays into a metastable state rather than the original state. In Pro.7Cao.3Mn03, 
a charge/orbital-ordered insulating state is melted and a ferromagnetic metallic state is induced 
by the photo-irradiation [5j[6]. Such a photo-induced insulator-to- metal transition has also been 
demonstrated in Ndo.sCao.sMnOs [7j. 

In this paper, to gain deep insight into the ultrafast non-equilibrium dynamics of the complex 
charge-spin-orbital state, we investigate the real-time dynamics of a charge-excited state in an 
e g -orbital model from the viewpoint of the time evolution of wavepackets by exploiting numerical 
techniques. We discuss the effects of the pair-hopping interaction on the propagation of holon 
and doublon. 
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Figure 1. Electron configuration of (a) ground state and (b) charge-excited state. Dotted lines 
indicate hopping connection between adjacent a. orbitals 

2. Model and numerical method 

Let us consider an e g -orbital degenerate Hubbard model on the linear chain along the z direction 
with N sites and one electron per site (quarter filling), described by 

Piaapipa' 

+ J Yj 4aa4/3v> d ic^'dil3* + J' 4^4^'^^, (1) 

where di Ta is an annihilation operator for an electron with spin a (=t>4-) i n orbital r (=a for 
3z 2 -r 2 ; p for x 2 —y 2 ) at site i, and Piro-=^J ro -^iTo-- The hopping amplitude is given by t aa =l 
(taken as energy unit), t a p=tp a =tpp=0. Note that due to the orbital anisotropy, 3z 2 — r 2 orbital 
becomes itinerant, while x 2 —y 2 orbital is localized. Regarding the on-site Coulomb interaction, 
U is the intra-orbital Coulomb interaction, U' the inter-orbital Coulomb interaction, J the inter- 
orbital exchange interaction, and J' the pair-hopping interaction between different orbitals. The 
relation U=U'+J+J' holds due to the rotation symmetry in the local orbital space and J=J' is 
assumed due to the reality of the orbital function [2], so that we have two independent interaction 
parameters. In this paper, we set U'=10 and study the dependence on J'. We use the unit such 
that h=l, and the time is measured in units of h/t aa . 

We investigate the dynamics of the charge-excited state from the viewpoint of the propagation 
of wavepackets [8]. For this purpose, we exploit density-matrix renormalization group (DMRG) 
techniques. First, the ground state \ipo) is obtained by an ordinary static DMRG method |9|. 
The finite-system algorithm is employed with the use of open boundary conditions. Then, as 
an initial state at time t=0, we prepare a charge-excited state by creating a holon-doublon pair, 
defined as 

i^) = E t -' c Lc i+ i T vi^o). (2) 

tt'ij 

The time evolution of \ifii d ) is computed by an adaptive time-dependent DMRG method |1(J|. II 1 j . 
The second-order Suzuki- Trotter decomposition is used to describe the time-evolution operator 
as the product of local time-evolution operators with a small time step At=0.02. We keep up 
to m=300 states in renormalization steps, and the truncation error is kept below 10~ 7 during 
the time evolution. 

3. Results 

In Fig. 1(a), we depict the electron configuration of the ground state. Since the electron hopping 
connects a orbitals in adjacent sites, a orbital is singly occupied in every site so as to gain 
kinetic energy. Localized 5=1/2 spins exhibit quasi-long-range antiferromagnetic order due to 
the antiferromagnetic exchange interaction. For the charge-excited state, a holon-doublon pair 
is created at the center of the chain, as shown in Fig. 1(b). We note that a orbital is doubly 
occupied in a doublon site. 
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(b) spin 
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(c) orbital 
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Figure 2. Time evolution of (a) charge, (b) spin, and (c) orbital densities at J7'=10 and J'=0. 
A holon-doublon pair is introduced at the center of the chain of iV=32. Time-series results from 
time t=0 to t=10 are arranged with shifted upward along the vertical axis. 



To gain insight into the dynamics of the charge-excited state, we measure the time evolution of 
the charge, spin, and orbital densities, defined by n c (i,t)=J2 T a(Pira)t, n s (i,t)=J2 T {Pirt-Piri)t, 
and n Q (i, t)=J2 a (Piacr—Pii3a)t, where (• • -)t denotes the expectation value using the wavefunction 
at time t. In Fig. 2, we present our numerical results at J'=0. Concerning the charge density, 
a holon-doublon pair is completely localized at the chain center at t=0, i.e., n c (N/2, 0)~2 and 
n c (N/2+l, 0)~0, as shown in Fig. 2(a). As the time evolves, holon and doublon wavepackets 
expand in the system. We see that the wavefront of the holon wavepacket propagates toward 
right, while that of the doublon wavepacket moves left, but there is no fine structure in the region 
between two wavefronts due to interference effects. Note that holon and doublon wavepackets 
propagate with the same speed, since holon and doublon move through the electron hopping of 
the same amplitude t aa . On the other hand, as shown in Fig. 2(b), the spin density is not affected 
by the charge excitation, since we assume Sf ot =0 for the charge-excited state, where Sf ot is the 
z component of total spin. As shown in Fig. 2(c), the orbital density is equivalent to the charge 
density. Once two electrons occupy a orbital in the doublon site, j3 orbital is always vacant due 
to the absence of the pair-hopping process at J'=0, so that n c (i,t)=n (i,t)=J2 a {Piaa)t- 

At finite J', however, a pair of electrons accommodated in a orbital can transfer to /3 orbital 
through the pair hopping, leading to a difference in the propagation of holon and doublon. In 
fact, we find that the doublon wavepacket propagates at half speed in comparison to the holon 
wavepacket, as shown in Fig. 3(a). Let us here discuss the effects of the pair hopping on the 
doublon propagation. Assuming that the electron pair occupies either a orbital or f3 orbital, 
the local two-electron state is described by the linear combination of |£aff) |ia,J.) and , 
where \iTa)=d iT(T \0) and |0) is the vacuum state. With these two states, the eigenstates are given 

by \ia)=(\ia\)\ial)—\iPt)\iPl))/V^ and \ib)=(\ic(t)\ial)+\i/3-f)\i/3l))/y/2 with its eigenenergy 
E a =U—J' and Eb=U+J', respectively. Among them, we suppose that the lowest-energy state 
\ia) represents the doublon state in the strong-coupling limit. Then, doublon and single electron 
can be exchanged by an effective hopping, 

(i + laa'\(ia\t aa (c\ a(7 Ci + i ar7 + h.c.)\iaa')\i + la) = --i a( A,-a'- (3) 



Namely, the hopping amplitude for doublon is reduced to half of the original hopping amplitude 




Figure 3. Time evolution of (a) charge, (b) spin, and (c) orbital densities at J7'=10 and J'=2. 



due to the pair hopping. On the other hand, holon and single electron are exchanged with the 
original hopping amplitude. Thus doublon propagates at half speed comparing with holon. 

For the spin state, no structure appears in the spin density even at J'=2, as shown in Fig. 3(b). 
Figure. 3(c) presents the orbital density. We can see a correspondence of the orbital density to 
the charge density, while they are not equivalent any longer. In fact, regarding the right-moving 
wavefront, the orbital density is quite similar to the charge density. The left-moving wavefront 
of the orbital density propagates with the same speed as that of the left-moving wavefront of 
the charge density. However, the difference of the orbital density from the uniform value varies 
between positive and negative as the time evolves, indicating that the electron pair actually 
occupies not only a orbital but also f3 orbital. Thus the orbital structure around doublon is 
dynamically deformed due to the pair hopping. 



4. Summary 

We have studied the time evolution of the charge-excited state of the one-dimensional e g -orbital 
degenerate Hubbard model by using DMRG techniques. We have shown that doublon propagates 
at half speed in comparison to the holon motion because of the pair-hopping process peculiar to 
multi-orbital systems. 
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